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In this work, we study the arising of correlations among some isoscalar ( K 0 , Q 0 , and I 0 ) and 
isovector (J, L 0 , K° yra , Q° ym , and 7° ym ) bulk parameters in nonrelativistic and relativistic hadronic 
mean-field models. For the former, we investigate correlations in Skyrme and Gogny parametriza- 
tions, as well as in the nonrelativistic (NR) limit of relativistic point-coupling models. We provide 
analytical correlations among bulk parameters for the NR limit, discussing the conditions in which 
they are linear ones. Based on a recent study [B. M. Santos et al, Phys. Rev. C 90, 035203 (2014)], 
we also show that some correlations presented in the NR limit are reproduced for relativistic models 
presenting cubic and quartic self-interactions in the scalar field a, mostly studied in this work in the 
context of the relativistic framework. We also discuss how the crossing points, observed in the den¬ 
sity dependence of some bulk parameters, can be seen as a signature of linear correlations between 
the specific bulk quantity presenting the crossing, and its immediately next order parameter. 
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I. INTRODUCTION 

In a general way, there are at least two differ¬ 
ent and competitive approaches in the treatment of 
nuclear matter. One of them is based on nucleon- 
nucleon interactions, from which many-nucleon mi¬ 
croscopic relativistic/nonrelativistic Brueckner-Hartree- 
Fock (BHF) [1] calculations are performed to obtain in¬ 
formation regarding the entire nuclear system. These 
calculations depend on the chosen nucleon-nucleon po¬ 
tential, classified as phenomenological (Reid, Urbana and 
Argonne interactions, for instance) or theoretical one- 
boson exchange ones (Paris, Bonn and Nijmegen inter¬ 
actions, for instance). For a short review, see Ref. [2]. 
Such nucleon-nucleon interactions reproduce experimen¬ 
tal data on phase-shifts and deuteron properties and are 
implemented in complicated many-nucleon BHF codes in 
order to obtain nuclear matter properties. Alternatively, 
a second (macroscopic) approach does not use nucleon- 
nucleon interaction itself. This is the case of nonrela¬ 
tivistic mean-field models like Skyrme [3] and Gogny [4] 
ones. In relativistic framework, on the other hand, the 
most used models are the relativistic mean-field (RMF) 
Walecka model [5], and its improved versions [6]. 

In both approaches, many of the models and approxi¬ 
mations used are intrinsically related to the fact that an 
analytical expression for the nucleon-nucleon potential is 
unknown. Therefore, in a theoretical point of view, cor¬ 
relations between two or more observables acquire enor¬ 
mous importance due to the fact that they reduce the 
set of independent relevant quantities to be used in the 
construction of nuclear models, avoiding redundant free- 
parameter fittings. In an experimental point of view, on 
the other hand, a correlation between two observables A 
and £>, for instance, allows the complete knowledge of B 


if A is experimentally well constrained. 

In the few-body nuclear physics, for example, the Tjon 
line [7], establishes a correlation between the binding en¬ 
ergies of 4 He and triton, B a and B t , respectively. The 
parametrization of the numerical results for several two- 
nucleon potentials [8] concludes that such a correlation 
reads roughly B a = 4.72 (B t — 2.48), in MeV units. It 
means that, if B t is calculated by using a two nucleon- 
nucleon potential, the value of B a is predicted, even be¬ 
fore any four-nucleon calculation. 

Concerning correlations among bulk parameters, there 
are still a few of them well established in the literature. 
One of them, usually known as the Coester line [9], cor¬ 
relates the saturation density p 0 and the nuclear mat¬ 
ter binding energy B a . It was analyzed in Ref. [9] in a 
two-nucleon model interaction, by varying its tensor force 
contribution while keeping the deuteron binding energy 
fixed. They showed that the function B 0 x p a roughly 
follows a line, in fact a band, also observing its similar¬ 
ity with the curve constructed from B 0 and p Q obtained 
from distinct two-nucleon potentials. Even modern cal¬ 
culations using Brueckner-Hartree-Fock, and including 
single-particle contribution in the continuum, change the 
results but preserve the Coester line. Another correlation 
was studied in Ref. [10], and involves the relationship 
between finite nuclei spin-orbit splittings and the ratio 
to* = M*/M for a class of finite-range (FR) RMF mod¬ 
els, being M* the nucleon Dirac effective mass at p = p 0 
and M the nucleon rest mass. The authors showed that 
such splittings are experimentally well reproduced if the 
range for m* is constrained to the following inequality, 

0.58 ^ to* ^ 0.64. (1) 

We will refer the above relation extracted from Ref. [10] 
as the FRS constraint. It is important because it re- 
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lates bulk infinite nuclear matter calculation with the 
closed shell finite nuclei energy spectrum. Still regard¬ 
ing relationships among quantities of finite nuclei and 
infinite nuclear matter, we also point out the correlation 
between the neutron skin thickness in 208 Pb, and the 
liquid-to-solid transition density in neutron-rich matter, 
investigated in Ref. [11] with RMF models. 


In a recent paper [12], we have provided analytical ex¬ 
pressions linearly correlating the symmetry energy (J), 
its slope (L 0 ) and curvature (K° ) at the saturation den¬ 
sity in the framework of the nonrelativistic (NR) limit of 
nonlinear point-coupling (NLPC) versions of the Boguta- 
Bodmer models [13] (FR-RMF models presenting cu¬ 
bic and quartic self-interactions in the scalar field a). 
From such analytical correlations, we were able to predict 
ranges for L a that Boguta-Bodmer models must present 
in order furnish good values for finite nuclei spin-orbit 
splittings. In this paper, we further investigate the aris¬ 
ing of correlations among bulk parameters of infinite nu¬ 
clear matter at zero temperature in nonrelativistic and 
relativistic hadronic mean-field models. For the former, 
we choose some Skyrme and Gogny parametrizations as 
well as the NR limit of NLPC models. For the latter, we 
choose FR-RMF models, specially the Boguta-Bodmer 
ones, and show the conditions to be satisfied in order to 
their bulk parameters present correlations. In particu¬ 
lar, we show, in a general way, how linear correlations 
are connected with crossing points in the density depen¬ 
dence of bulk parameters. For this purpose, we generalize 
the procedure used in Ref. [14], where the authors associ¬ 
ated the crossing point in the incompressibility function 
for different nonrelativistic Skyrme models, but not for 
FR-RMF ones, with the linear correlation between K a 
and Q 0 (incompressibility and skewness coefficient at p a , 
respectively). We show here that there are other cross¬ 
ing densities different from p c ~ 0.7p o , value found in 
Refs. [14, 15], for different bulk parameters in nonrel¬ 
ativistic models, as well as in relativistic ones. In the 
following, we will see that such crossing densities may be 
seen as signatures of linear correlations between higher 
order derivatives of the specific bulk parameter present¬ 
ing the crossing. 


The paper is organized as follows. In Sec. II we show 
how linear correlations are connected with crossing den¬ 
sities in bulk parameters of infinite nuclear matter. In 
Sec. Ill, we use some results of our previous study of 
the NR limit [12], in order to apply the calculations of 
the previous section. In Sec. IV, we use the predictions 
of the NR limit to study the conditions that establish 
linear correlations in FR-RMF models. In Sec. V, we 
study relativistic and nonrelativistic models presenting 
two isovector coupling constants in the context of the 
correlation between the symmetry energy and its slope. 
Finally, we present our main conclusions in Sec. VI. 


II. CONNECTION BETWEEN LINEAR 
CORRELATIONS AND CROSSING DENSITIES 


In the literature, it is verified that crossing points can 
occur in the density dependence of the symmetry en¬ 
ergy [16], in pure neutron matter equation of state [17], 
and in pairing gap of nuclear matter [18], for instance. 
Recently, in Refs. [14, 15], the authors found a specific 
crossing point in the incompressibility of nuclear matter 
for different nonrelativistic Skyrme models, and not con¬ 
firmed in FR-RMF ones. They showed that this cross¬ 
ing density (p c ~ 0.7p o ~ 0.11 fm -3 ) when used in 
the calculation of the derivative of the incompressibil¬ 
ity, makes this quantity better correlated to the centroid 
energy of the isoscalar giant monopole resonance than 
K a . The authors also pointed out that this crossing den¬ 
sity p c is closer to the average density in the 208 Pb nu¬ 
cleus, (p) = 0.12 fm~ 3 , than the saturation density it¬ 
self, p 0 ~ 0.16 fm -3 , with (p) also obtained by using 
the Skyrme model. The question we pose here is what 
properly mean, or indicate, such crossing points. In the 
following, we try to answer this question. 

Actually, crossing points in different bulk parameters 
of infinity nuclear matter can be viewed as a signature 
of linear correlations between higher order derivatives of 
that particular bulk parameter presenting the crossing in 
its density dependence. In order to make it clear, we pro¬ 
ceed to generalize the calculation performed in Ref. [14], 
where the authors associated a crossing in the K ( p ) func¬ 
tion with the linear correlation between K 0 and Q 0 for 
some Skyrme parametrizations. Firstly, let us define a 
function of the density, J-'(p), expanded in terms of the 
dimensionless variable x = p Z Po , and around the satura- 
tion density as, 


^ ^ ^ x , F"{Po) 2 , r"( Po ) 3 

F(p) = J-(p 0 )+J r ( Po)X-\ - ——x -I-—- x 
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with the derivatives of T{p) given by 
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(3 p) 


d_ r T^-^jpy 
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for m = 1 , 2 ,3 ... 

For the values m = 1,2, 3, for instance, and by noting 
that — = 3ir + 1 and -M -, we see that 
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and 


T"'{p) = (3 p) 
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F'(P) 


dp [ (3p ) 2 _ 
(3p ) 3 d 3 T 


= (3 pf 
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d 2 T 


(3 p 0 ) 3 dx 3 
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= (3 ^ + 1 ) 3 ^T- ( 6 ) 


The pattern verified in Eqs. (2)-(6) allows us to write 
_F(p) and its derivatives in a compact form as, 


= (3z + l) m {^ m \po) + F {rn+1 Xpo) 

m+2 \p 0 ) 
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,( 7 ) 


for to = 0 , 1 , 2 ,3... and with T^ m ^{p 0 ), E^ m+1 \p 0 ), 
J r ( m+ 2 )( ( o 0 )... being the bulk parameters evaluated at the 
saturation density. 

In order to make our analysis simpler, we consider 
that /{3x + l) m is well described by its expansion 
until order x 3 . Now, let us assume that the bulk pa¬ 
rameters of the higher order derivatives of , namely, 
jr(m+i)( (0o ), E^ m+2 \p 0 ) and J r( - m+3 \p 0 ), are correlated 
with T^ipo), more specifically in a linear way, as given 
below, 

^ rn+1 \p 0 )=b 1 + a 1 ^ m \p 0 ), ( 8 ) 

E^ m+2 \p 0 ) = b 2 + a 2 ^ m \p 0 ), (9) 

^ m+3 \p 0 ) = b 3 +a 3 ^ m \p 0 ), ( 10 ) 

with cii and bi independent of any other parameters of 
the model. If such linear correlations hold, then, Eq. (7) 
can be rewritten as 


E {m \p) ~ (3a; + l) m x 

x | f(x)^ m) (p 0 ) + bix + ^x 2 + ^a; 3 |, (11) 

where 

f(x) = l + aix + ^x 2 + ^x 3 . ( 12 ) 

For a specific value (or values) of the density, named as 
p c , that makes f(x c ) = 0 , with x c = pc 3p po , E^ m \p c ) will 
be exactly the same for any hadronic model, according to 
Eq. (11) since the bi parameters are model independent. 
Therefore, the functions T^ m \p) of any model will cross 
each other exactly at p = p c . This is the crossing point 
presented in the (p) function. Thus, one can see this 
crossing as a signature of the linear correlations given 
by Eqs. (8)-(10). Naturally, such a signature only holds 
if we can express the function E^ m \p) in terms of the 
expansion presented in Eq. (7). 

In the case of jd m )/( 3 a; -f l) m is expanded until or¬ 
der x 3 , Eq. (12) generates a cubic equation, when one 
imposes f(x c ) = 0 , in order to localize the values of 
p c . In the general case, in which is needed to expand 
J"( m )/(3x + l) m until order x N , then 

f( Xc ) = i+E^r = °’ ( 13 ) 


produces an equation of order N to be solved in order to 
determine the possible values of x c , and, consequently, 
the crossing densities p c of the function E^ m \p). 

The determination of linear correlations from the 
searching of crossing points is suitable, for instance, for 
the isovector bulk properties of hadronic models, in which 
the quantities and its derivatives can be expanded around 
the saturation density exactly as in Eq. (7). The en¬ 
ergy per particle of a system of proton fraction given by 
y = p p /p ( p p is the proton density) can be expanded in 
terms of the isospin asymmetry parameter /l = 1 — 2y as 

E(p, ft = E{p) + S(p)(3 2 + S 4 (p)/3 4 + ■■■ (14) 

where E(p) is the energy per particle related to the sym¬ 
metric nuclear matter (/3 = 0). The other two coeffi¬ 
cients of the expansion are, respectively, the symmetry 
energy and the fourth-order symmetry energy, with the 
density dependence, associated with the isovector sector, 
expanded in terms of the density as 

S(p) = J + L 0 x+^x 2 + ^x 3 + --- (15) 

and 

S 4 (p) = J 4 + L\x + ^± x 2 + + ... ( i 6 ) 

The derivatives of S(p) are defined exactly as in Eq. (3), 
i. e., L = 3pf|, K sym = (3p) 2 0, Q sym = (3p) 3 0, and 
so on. The isovector bulk parameters, L 0 , K° ym and Q° ym 
are the derivatives evaluated at p = p 0 , and J is given 
by J = S(p 0 ). Analogous quantities are also defined 
from S 4 (p). Based on this structure, the search of linear 
correlations from the location of crossing densities can be 
naturally performed. This will be done for nonrelativistic 
and relativistic models in the next two sections. 

As a last remark, we mention here that the procedure 
described above was first used in Ref. [14] specifically to 
justify the crossing point in the K(p) x p curve, and not 
as a route to find linear correlations as we are doing in 
the present work. Moreover, we are also generalizing this 
method to any bulk parameter. 


III. CORRELATIONS IN NONRELATIVISTIC 
MODELS 


A. Theoretical framework of the NR limit 

In the background of nonrelativistic mean-field models, 
we analyse some parametrizations of Skyrme and Gogny 
models, and also those from the NR limit of NLPC ver¬ 
sions of the Boguta-Bodmer model. It is important to 
mention that as the FR-RMF models, relativistic point¬ 
coupling ones also describe very well the infinite nu¬ 
clear matter bulk parameters and finite nuclei proper¬ 
ties [19-24]. In Ref. [24], for instance, the authors were 
able to obtain, by using a NLPC model, ground state 
binding energies, spin-orbit splittings, and rms charge 
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radii of a large set of closed shell nuclei, as well as, of 
nuclei outside the valley of beta stability, clearly showing 
the success of these kind of model. Their nonrelativistic 
versions, besides follows this same pattern, at least con¬ 
cerning the infinite nuclear matter that is the scope of 
our work, are also useful in the sense that they can be 
used to predict correlations also exhibited in FR-RMF 
models, as observed in our previous study of Ref. [12]. 

The relativistic NLPC versions of the Boguta-Bodmer 
models are described by the following Lagrangian density 


The symmetry energy iS (NR) is used in order to obtain its 
slope, curvature and skewness. The results are 

e™ if) = ^ (m + \gI„ +sv +1 v) 

+ 3 G 2 v p, (23) 

^ i~ M + 5 Glp + 40 V + 132 Bp 3 ) , (24) 


£nlpc = - M)ijj - + ^Gs(V’V') 2 

+ + j(#) 4 - ^GtvWV^) 2 , (17) 


that mimics the two-, three- and four-body pointlike in¬ 
teractions. In this equation, the last term is included 
in order to take into account the asymmetry of the sys¬ 
tem (different number of protons and neutrons). In the 
NR limit of the NLPC model, and by using the mean- 
field approximation, the energy density functional at zero 
temperature for asymmetric nuclear matter is written as 


e^(p,y) = (G 2 v ~G 2 s )p 2 ~Ap 3 -Bp 
+ G 2 v p 2 (2y — l) 2 + 


Ap3,(18) 


10 M*(p,y) 


where the effective mass is 
M*{p,y) = 


M 2 


(M + G%p + 2 Ap 2 + 3Bp 3 )Hs_ 


(19) 


and 

QZ\p) = (m \g 2 s P + 20V + 165v) , 

(25) 

respectively. 

The coupling constants of the model are G|, G 2 , A, B 
and G 2 v . The first four of them are adjusted in order to 
fix p 0 , B 0: K 0 and M*. This is done by solving a sys¬ 
tem of four equations, namely, e (NR) (p 0 ,1/2 )/p 0 = — B 0l 
IF <NR) (p 0 ,1/2) = K a , P (NR) (p 0 ,1/2) = 0 (nuclear mat¬ 
ter saturation), and M*(p 0 ,1/2) = M*. The last cou¬ 
pling constant, G 2 V , is obtained by imposing upon the 
model the requirement of presenting a particular value 
for J = £> (NR) (p 0 ). The explicit forms of the constants 
G 2 , G 2 , A, B and G 2 V in terms of to*, p 0l B a , K a and 
J, are given in the Appendix. 


B. Results from the isovector sector 


with Hh = 2a [ys + (1 — y)l], and A = (37r 2 /2)3, see 

Ref. [12]. 

From the energy density in Eq. (18), it is possible to 
obtain pressure, incompressibility and the symmetry en¬ 
ergy of the model, since P = p 2 a( ' £ Q p p ' ) , K = and 


c _ l [ d 2 (e/p) 
° ““ 8 dy 2 


These expressions are, respectively, 


given by 


P mK) (p,y) = (G 2 - G 2 )p 2 - 2 Ap 3 - 3 Bp 4 
+ G 2 v p 2 (2y-l) 2 + ^x 

x ^ M + -G 2 p + 8 Ap 2 + —-Bp 3 ^ 


pi, (20) 


AT (NR) (p, y) = 18(G 2 - G 2 )p - 54Ap 2 - 108Bp 3 
+ 18G 2 v p(2y - l) 2 + x 

x + 4G 2 p + — Ap 2 -\ — —Bp 3 ^j p 3 , (21) 


and 


Aps 

6Af*(p, 1/2)' 


( 22 ) 


In Ref. [12], we rewritten the coupling constants of the 
model in terms of the bulk parameters to*, p Q , B 0 , and 
K 0 (an analogous procedure is done in the context of the 
Skyrme models in Ref. [25]). This method allowed us to 
explicitly write the slope of the symmetry energy at the 
saturation density also as a function of m*, p D , B 0 , and 
K 0 , and thus, find the following correlation between J 
and L 0 , namely, 

L 0 = 3J + b(m*,p 0 ,B 0 ,K 0 ), (26) 


with 


b(m*,p 0 ,B 0 ,K 0 ) = 

x(±^(3M 2 - 
1 9 to* 


1 

(3 M 2 - 19 E°M + 18 Ef) * 

1AM E°) - 5 E°(M 2 - 5ME°) 

+ 30 B 0 Ef --^E^M\, (27) 


where E° = 3A pi /10 M. 

The expression in Eq. (24) along with the correlation 
in Eq. (26), and the definitions of the coupling constants 
in terms of the bulk parameters are used to find 

CU = (Lo - 3 J)p(po) + q(Po , B 01 K 0 ) (28) 


S< NR >(p) = G 2 v p + 
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where 


P(Po) 


3 (5 M 2 - 18 E°M) 
3 M 2 - UME° 


(29) 


and 


q(Po, B 0 , K 0 ) 


5 M 2 E° + 90 B 0 ME° - ?§-K 0 E°M 
3 M 2 - UME° 


,(30) 


As pointed out in Ref. [12], if J and K„ are kept fixed 
in Eq. (28), K° yin will present a linear correlation with 
L„, since, the binding energy and the saturation den¬ 
sity are well established closely around the values of 
B 0 = 16 MeV and p 0 = 0.15 fm -3 , and do not vary 
too much for each parametrizations of any nuclear mean- 
field model. For this reason, we consider B 0 and p a as 
constants hereafter. The same conditions also occur with 
the skewness of S of the NR limit, due to its analytical 
structure in Eq. (25). This quantity, obtained through 
Q ° rm = Q™(p 0 ), can a l so be written in the following 
form, 


Q° ym = (Lo - 3 J)u(p 0 ) + v(p 0 , B 01 K a ) 

with 

f ^ 5 (3M 2 - 22 E°M - 144 E° 2 ) 

U ^ Po ’ ~ 3 M 2 - 14 ME° 

and 


(31) 


(32) 


v(po,B 0 ,Ko) — 


1 


3 M 2 - UMES,' 


850 


x { 5 E°M (. M + 52 E°) - —K 0 E°M 

+ 150 B 0 E° (9 M - 8 E°) !>. (33) 


The linear dependence of K° ym and Q° ym as a function 
of L a for fixed J and K a is depicted in Fig. 1. 

We remark to the reader that it was possible to in¬ 
vestigate how K° and Q° ym depend on L a , by using 
parametrizations presenting J fixed and different val¬ 
ues of the effective mass, specifically in the range of 
0.50 ^ m* ^ 0.80. In this way, it was possible to keep 
J fixed and still be able to vary L a from the variation of 
b(m*,p 0 , B a , K 0 ) in Eq. (26). 

The increasing of K° yra , and decreasing of Q° yrn as a 
function of L a , are also verified for Gogny interactions 
presented in Ref. [26]. In this work, the authors anal¬ 
ysed the isovector properties of this specific nonrelativis- 
tic model, providing analytical expressions for symmetric 
and asymmetric nuclear matter, see Fig. 2. Notice that 
although some parametrizations present a linear depen¬ 
dence, this behavior is not verified for all of them. 

According to the discussion of Sec. II, the searching of 
such linear correlations could also have been done if we 
had looked for their possible signature, in this case in the 
density dependence of the symmetry energy slope. If the 



FIG. 1. (a) K° yra and (b) Q° ym as a function of L 0 for the NR 
limit for different values of J. The effective mass varies in the 
range of 0.50 ^ m* ^ 0.80. For each panel, p 0 = 0.15 fm -3 , 
B 0 = 16 MeV and K 0 — 270 MeV. 
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FIG. 2. (Color online) (a) K° yra and (b) Q° ym as a function 
of L 0 for Gogny parametrizations studied in Ref. [26]. 


function L (NR) (p) can be expanded, and if its density de¬ 
pendence presents a crossing point, then one can ensure 
at least the linear behavior presented in Fig. la. In fact, 
this crossing around p^./po = 0.47 is verified in Fig. 3 in 
the NR limit for parametrizations in which the values of 
J and K 0 are kept fixed. 

It is worthwhile to note that the crossing density dis¬ 
played in Fig. 3, namely, p^/po = 0.47, is exactly the 
same for all curves. In this figure, we group three 
different sets of parametrizations, each one presenting 
three distinct values for the pair ( J , K 0 ). The val¬ 
ues were chosen inside the ranges of 25 ^ J ^ 35 MeV, 
and 250 ^ K a ^ 315 MeV. One can still notice that 
the quantity L (NR> (p R ) is different for each set. How¬ 
ever, the values of L (NR) (p R ) of these parametrizations 
present an overlap of 48% with the constraint estab¬ 
lished in Ref. [27], namely, T(p R ) = 47.3 ± 7.8 MeV. We 
still point out that the range of 25 ^ J ^ 35 MeV is 
showed to be totally compatible [6] with experimental 
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FIG. 3. (Color online) Symmetry energy slope as a function of 
p/po for the NR limit, Eq. (23) . In this figure, p a = 0.15 fm -3 , 
and B 0 = 16 MeV. 


values from analyses of different terrestrial nuclear ex¬ 
periments and astrophysical observations [28], and the 
range of 250 ^ K a ^ 315 MeV was based on the recent 
reanalysis of data on isoscalar giant monopole resonance 
energies [29]. 

If we consider the expansion of L (NR) /(3a; + 1) until 
order x, namely, L(p) ~ (3x + 1 )(L a + K° ym x), it is 
possible to explain the crossing point in Fig. 3 if K° yxa 
is linearly correlated with L 0l i. e., if we can write 
^C°ym = bi + a\L 0 , and if there exist a real root for 
the equation f(x c ) = a\x c + 1 = 0, since in this case 
one writes L(p) = (3a: + 1) {L a (a\X + 1) + bix}. As 
we can see in Eq. (28) and in Fig. la, the linear cor¬ 
relation between K° and L 0 holds in fact for J and 
K 0 fixed. However, if we solve the equation /( x c ) = 0 
we find /9 R /po = 0.41 (since Gq = p(p 0 ) = 5.13 and 
Pc/Po = 3x R + 1), not exactly the same crossing den¬ 
sity presented in Fig. 3. This suggests that /( x c ) should 
be modified in order to produces a more exact root, what 
means that L (NR) /(3x + l) actually needs to be expanded 
in higher orders in x , at least in the density region around 
the crossing density p R /po = 0.47. If we now take the 
expansion of L (NR) /(3x + l) until order x 2 , the crossing 
point will exist if K° yla and Q° ytn are linearly correlated 
with L 0 . The latter correlation is verified in Eq. (31) 
and Fig. lb as we already discussed. For this case we 
will have 


L{p) ~ (3x + 1) \ L 0 + K° x + 


Q° ym 2 


(3x + 1) | L 0 {^-x 2 + aix + i) + Y x + 



with bi = q — 3Jp , cq = p, 62 = v — 3Ju , and 02 = u. 
The crossing is explained if /(x R ) = 1 + px R + f x r2 = 0, 


15 M 2 - ME°M- 


c 5 (3M 2 - 22 E°M - 144 Ef) 

(315M 4 - 2700F;“M 3 +1676 E° 2 M 2 + 20160 EfM) 2 \, 


(35) 


that produces p R /po = 0.46, a value much more close to 
the crossing density than p R /po = 0.41, found previously. 

Therefore, it becomes clear that a crossing point in 
a density dependence of a bulk parameters indicates a 
route for the searching of linear correlations in its higher 
order derivatives. Nevertheless, we point out to the 
reader that a crossing point itself does not ensure lin¬ 
ear correlations in all higher order bulk parameters. For 
the previous analysis, for example, one can not affirm 

that L a will be correlated with 1° = (3p 0 ) 4 () , 

\ p J p= Po 

simply by the fact that the expansion of i (NR) /(3x + 1) 
until order x 3 is better than those until order x 2 . It is 
needed to check whether such linear correlation really 
holds. Actually, the crossing ensure at least the linear 
correlation between L 0 and the immediately next order 
bulk parameter K° ym . As a matter of fact, we investigate 
if L 0 correlates with 1° m by obtaining the expression of 
the fourth order derivative of <S (NR) , namely, 


I^\p) 


28A pi 
3 M 2 


(^—M + -G 2 p 


20 . 2 330 „ 

—Ap 2 + — Bp 



(36) 


From Eq. (36), is possible to write I° ytn = 7 ( ” n R) (p c) ) in 
terms of L 0 , J, p D , B 0 and K a in the following way, 


Cn = (Lo - 3 J)r{ Po ) + s( Po , B 0 , K 0 ) 

with 

^ a 20 (3 M 2 - 56 ME° + 324 Ef) 

r ( Po > = 3 M 2 - 14ME° 

and 


(37) 


(38) 


s(po,B 0 ,Ko) — 


3 M 2 - UME°' 


1080 B 0 E° (M - 2 E°) 


- 12 E°M (M - 18 E°) - ^-K a E°M (39) 


Thus, it is verified that I s ° m also linearly correlates with 
L a , under the same conditions that makes L a also corre¬ 
lated with K° ym and Q° yra , namely, fixed values for J and 
K 0 . Furthermore, a new expansion of L (NR) /(3x + 1) 
until order x 3 generates the cubic equation /(x R ) = 
1 + px\ + % x c 2 + ii^c 3 = 0 , presenting a root corre¬ 
sponding to Pc/Po = 0.47 (since r(0.15) = 14.17), the 
exact value for the crossing density. 
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Following these same ideas, we search for signatures of 
linear correlations in the Skyrme model. At this point, 
we remind the reader there is no unique crossing point at 
the density dependence of the symmetry energy and its 
slope for the Skyrme model, as we can see in Figs. 4 and 5, 
respectively, where we display the 240 parametrizations 
of Ref. [3] . 



FIG. 4. (Color online) Symmetry energy as a function of p/p 0 
for the 240 Skyrme parametrizations of Ref. [3] . 



FIG. 5. (Color online) Symmetry energy slope as a function 
of p/po for the 240 Skyrme parametrizations of Ref. [3]. 

This lack of a unique crossing in the density dependence 
of S(p) and L{p) functions can also be seen in Fig. 2 
(left) of Ref. [30], where the authors studied 21 Skyrme 
parametrizations. 

Specifically for the density dependence of the symme¬ 
try energy slope, we found a crossing density for the 
SV [31], SkO [32], SkO’ [32], SkT3 [33], SkRA [34] and 
Ska35sl5 [3] parametrizations at p/./ po = 0.38. It is dis¬ 
played in Fig. 6. 

According to the discussed so far, these Skyrme 
parametrizations will present linear correlation at least 
regarding K° and L 0 . This is confirmed in Fig. 7a. 
Moreover, in Figs. 7b and 7c it is also checked the 



FIG. 6. (Color online) Symmetry energy slope as a function 
of p!po for some Skyrme parametrizations. 


linear correlations of L a with Q° ym and J s ° m , respec¬ 
tively, like in the case of the NR limit. For the sake 



FIG. 7. (Color online) (a) K° ym (b) Q° ym , and (c) J° ym as a 
function of L 0 for the Skyrme parametrizations of Fig. 6. The 
linear fits are indicated in each panel. 


of completeness, we have checked that the expansion 
of L(p)/(3x + 1) that approaches to the exact function 
around p/p 0 — 0.38 is taken until order x 3 for these 
Skyrme parametrizations. Therefore, the angular coef¬ 
ficients found in Fig. 7 are used to define the function 
f(Xc) that needs to be null. This will provide the cubic 
equation 1 + 3.76a;]; — = 0, that has one 

of the roots given by pl/po = 0.38, precisely the crossing 
density verified in Fig. 6. 

Still at the framework of the Skyrme parametrizations, 
another crossing point is observed in the isovector sec¬ 
tor, specifically in the density dependence of the sym¬ 
metry energy itself. As pointed out in Fig. 8, such 
a crossing occurs at pi/po = 0.89. We have noticed 
that for these parametrizations, J is linearly correlated 
with L ni K° , Q° and 1° as one can see in Figs. 9 

u ' sym" * sym sym o 

and 10. The angular coefficients of these lines are used 
to define the quartic equation to be solved in order 














P'Po 


FIG. 8. (Color online) Symmetry energy as a function of p/po 
for some Skyrme parametrizations. 


to determine the value of the crossing density, namely, 
1 + 28.8a:® + if ^xf - if^xf + ifix® 4 = 0. One 
root of this equation provides the value pi/po = 0.89, 
observed in Fig. 8. 



FIG. 9. (Color online) (a) L 0 and (b) 77° ym as a function of 
J for the Skyrme parametrizations of Fig. 8. The linear fits 
are indicated in each panel. 


As a remark, it is worthwhile to note that a linear 
correlation itself between two bulk parameters is not a 
sufficient condition to guarantee a crossing point in the 
density dependence of the immediately preceding bulk 
parameter. As an example of this statement, we focus 
on the analytical structure of the NR limit to find other 
two specific linear correlations. From Eqs. (28) and (31) 
it is straightforward to obtain 

Q° ym = -K ym + v- —■ (40) 

p p 

From Eqs. (28) and (37), an analogous expression can be 
written relating 7° and 7f 3 ° ym , namely, 


7° 


r 

= -K° 


rq 

P ' 


(41) 



FIG. 10. (Color online) (a) Q° ym and (b) 7° ym as a function of 
J for the Skyrme parametrizations of Fig. 8. The linear fits 
are indicated in each panel. 


Therefore, one see for a fixed value of K 0l that the bulk 
parameter 77 3ym is linearly correlated with Q/ ym and 7 3ym , 
see Fig. 11. 



FIG. 11. (a) Q°ym and (b) 7° ym as a function of A° ym for 
the NR limit. The effective mass varies in the range of 0.50 ^ 
m* ^ 0.80. For each panel, p 0 = 0.15 fm~ 3 and B 0 = 16 MeV. 

Here, we highlight that one can vary K° m without any 
information regarding L a . Let us remind that K° can 
also be written only as a function of the isoscalar param¬ 
eters, according to Eq. (12) of Ref. [12]. Thus, for fixed 
values of K 0 , it is possible to investigate how Q° ym and 
I“m depend on K° ynl only from the variation of m*. 

Notice that, the expansion of 7f 3 ( yl ^ ) /(3x + l) 2 until 
order x 2 can describe the exact function in Eq. (24) in 
a density region of subsaturation densities, see Fig. 12a. 
In spite of this, the quadratic equation constructed from 
the angular coefficients extracted from Fig. 11, namely, 
1 — 0.88xc sym + pXc sym2 = 0, presents no real roots, 
indicating no crossing points in the TG™* ip) function. 
This finding is confirmed in Fig. 12b. 
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FIG. 12. (Color online) as a function of p/p 0 - (a) Exact 

function for m* = 0.6, Eq. (24), compared with its expansion, 
(b) Exact function for different effective mass values. The 
bulk parameters fixed in the two panels are K 0 = 270 MeV, 
po = 0.15 fm -3 and B 0 = 16 MeV. 


Another example of linear correlations between bulk 
parameters and the lack of crossing points in the density 
dependence in one of them, is the case of the quanti¬ 
ties L 0 and J. According to Eq. (26), a linear corre¬ 
lation between J and L a is established if the function 
b(m *, po, B a , K 0 ) is kept fixed, i. e., if the analyzed 
parametrizations have exactly the same isoscalar bulk pa¬ 
rameters. Thus, parametrizations with different J values 
but same isoscalar parameters, present the behavior de¬ 
picted in Fig. 13a. 



FIG. 13. (a) L 0 x J correlation obtained from Eq. (26). (b) 
density dependence of symmetry energy, Eq. (22). 


studied can be approximated by its expansion, and si¬ 
multaneously, if the equation f(x c ) = 0 present nonzero 
real roots in the analysed range of densities. In the case 
of the NR limit, the latter condition is not satisfied in the 
study of the Kf R \p) function, implying in no crossing 
points in its density dependence in the range of subsatu¬ 
ration densities. 


C. Results from the isoscalar sector 


Regarding the quantities related to the isoscalar sec¬ 
tor of the nonrelativistic hadronic models, we underline 
here the relationship between K 0l Q a and I 0 . For the 
former two quantities, a correlation was firstly found in 
Refs. [14, 15] for some Skyrme parametrizations. In par¬ 
ticular, they found it as a linear one. Here we proceed 
to find, in the NR limit framework, the conditions that 
the parametrizations must satisfy in order to gives rise 
to the same relationship. For this purpose, we first need 
to obtain Q(p, y) from the energy density, Eq. (18). This 

r g 3 (g/p) ~ 

dp 3 


is done by calculating Q = (3 p) 3 
expression given by 


with the full 


Q (NR) (p,y) = —162Bp 3 + 


12A H, 


5/3 


5 M 2 


P 3 X 


x | M- -G 2 s p + 20 Ap 2 + 165 Bp 3 ) , (42) 

where Q 0 is defined by Q 0 = Q (NR) (p 0 ,1/2). 

By using the coupling constants in terms of the bulk 
parameters in the expressions of K a and Q 0 , it is possible 
to find the following relationship between these quanti¬ 
ties, 


where 


Qo — 'i(Po)K a T j {poi , nr ), 

3 (9 M 2 - 73 E°M + 90 Ef) 


i(Po) = 


and 

j{p 0 ,B 0l m*) = 


3 M 2 - 19 E°M + 18 Ef 


3 M 2 - 19 E°M + 18.Ep 2 


(43) 

(44) 


24 E\ 


^ (M 2 - 40ME°/3 + 60F;° 2 ) 


- 162R 0 (3 M 2 - 25 ME° + 40 E° 2 ) 


As one can see in Fig. 13b, these parametrizations do not 
generate any crossing point in the density dependence of 
S(p). From the point of view of the discussion of Sec. II, 
one can understand the lack of crossing points from the 
equation 1 + 3cc 3 = 0 , satisfied only for x 3 = —i. e., for 
pl/po = 0 (since we have p/po = 3x + 1). 

In summary, one can associate linear correlations as 
signatures of crossing points only if the exact function 


- 108F1° (M 2 - 10 ME° + 27 Ef) j. (45) 

From the structure presented in Eq. (43), we notice that 
the effective mass plays a crucial role for the correlation 
between Q 0 and K 0 . For fixed values of m*, this correla¬ 
tion is a linear one. We remember the reader that p a and 
B 0 vary only in a very narrow range around 0.15 fm —3 
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and 16 MeV, respectively. Thus, these bulk parameters 
can be considered constants for the hadronic mean-field 
models. For different to* values, Eq. (43) generates par¬ 
allel lines as indicated in Fig. 14a. The same linear corre¬ 
lation is also observed in some Gogny parametrizations, 
as indicated in Fig. 14b. 


bulk parameter associated to the immediately next order 
derivative of K(p). Nevertheless, for the NR limit, we 
can also check analytically if the next bulk parameter, 
J D , correlates with K 0 . From Eq. (18) and the definition 


I = (3pf 


9 J (£/p) 

dp 3 


, one obtains, 


O' 



'■•V 
/ •/ 

d > — m* = 0.5 (/ = - 2878) 

-- m* = 0.6 (j = -2926)' 
m* = 0.7 O'= -2961) 

- ■ ■ m* = 0.8 0 = - 2986)' 

. 

255 270 285 300 315 200 220 240 260 280 

K (MeV) K (MeV) 


' 1 ' 1 ' 1 

(b) 

i T 

V 

J_i_ 

/ - 


- 

O.V / 


/ • 

D1 

m 

DIM 

/ ♦ 

D1S - 

/ ▲ 

D260 

-/ ◄ 

_._i_i_i_i_i_ 

D280 “ 

_1_i_ 


-250 

-300 

-350 

-1-400 
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FIG. 14. (Color online) Q 0 as a function of K 0 for (a) 
the NR limit with different m* values, and for (b) some 
Gogny parametrizations studied in Ref. [26]. In panel (a), 
po = 0.15 frrR 3 and B 0 = 16 MeV. 

From the perspective addressed in Sec. II, the linear 
correlation between K 0 and Q 0 could also be sought, by 
searching for a possible crossing point in the density de¬ 
pendence of the incompressibility. In fact, as pointed out 
in Fig. 15, there are two of them, at p* / po = 0.21 (not 
shown), and p*/po = 0.79. Furthermore, notice that the 


r(NR) i \ 84A77 5 /3 2 

1 = 


5^2 20 2 330 3 

-M + -G s p —“I— —Bp' 


(46) 


From this expression, the fourth order derivative of the 
energy per particle evaluated at the saturation density, 
I 0 = / (NR) (p 0 ,1/2), can be written in terms of p a , B a , 
to*, and K 0 as 


I 0 = w{p 0 )K 0 + z(p 0 , B a , m*), 

where 

20 (25 E°M - 54 E° 2 ) 
W ^ Po ’ ~ ~ 3 M 2 - 19 E°M + 18 E° 2 ’ 

and 


(47) 


(48) 


z (poi Bo, 


X 


m ^ 3 M 2 - 19 E°M + 18 E° 2 X 

I 40E£ , 2 _ 56ME o + 3 24E?) 

TO* 

+ 1080R o (9 ME; - 32E° 2 ) 

- 72 E° (AM 2 - 77 ME° + 324 E° 2 ) 


(49) 



FIG. 15. (Color online) Incompressibility as a function of 
p/po for the NR limit, Eq. (21). 


the second point is quite close to the crossing density 
Pc/Po = 0-71 observed also for the Skyrme parametriza¬ 
tions studied in Ref. [14], 

Since crossings in the iC <NR) (p) function were found, 
they indicate a signature of a linear correlation, in this 
case at least between K a and Q 0 , since the latter is the 


Notice that once more, the effective mass needs to be con¬ 
stant for the parametrizations in order to ensure a linear 
dependence between I Q and K a , with angular coefficient 
given by w(0.15) = —4.16. 

For the sake of completeness, we use the angular coef¬ 
ficients i(p 0 ) and w{p 0 ) to calculate the crossing density 
in Fig. 15. First, we consider the energy per particle of 
symmetric nuclear matter as 

E(p) ~E 0 + ^ x 2 + Q?-x 3 + ^jz 4 , (50) 

then, the corresponding expansion for K ( p ) reads 


K ( p ) = 18p^ 
dp 


9 p z 


,d 2 E 


dp 2 

= 6(3cc + l)-w- + (3a; + l ) 2 


dE 


dx 


d 2 E 

Ih? 


~ (3x + 1) x 

K 0 + ( 9 K a + Qo)x + UQo + ~ ) a ; 2 -I— —3 


5 1 0 


(51) 

.(52) 


This expansion is observed to be consistent with the exact 
function, Eq. (21), as showed in Fig. 16. 
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p/Po 

FIG. 16. (Color online) Density dependence of the incom¬ 
pressibility compared with its expansion for the NR limit. 


It is worth noting that despite the extra term in 
Eq. (51) compared with T" in Eq. (5), the final ex¬ 
pansion, Eq. (52), is analogous to the general function 
F^ m \p) in Eq. (7). Therefore, all the procedure devel¬ 
oped in Sec. II also applies in the analysis of correla¬ 
tions and crossing points for the AT (NR> (p) function in the 
isoscalar sector. Indeed, this was done for some Skyrme 
parametrizations in Ref. [14]. From this point of view, it 
is possible to use Eqs. (43) and (47) to rewrite Eq. (52) 
as 


K(p) = (3x + l)x 



9 5 w , 

( 6l+ 2) 

* 2 + -yz 3 


+ j + (pj + |) X + ^-x 2 a; j . 


(53) 


Thus, the crossings points in Fig. 15a are justified if 1 + 
(^+i)x l ^-\-{ l 6i+w/2)xp 2 + (bw/2)x^ 3 =0. Two solutions 
of this cubic equation are p^/po = 0 . 21 , and p^/po = 

0.79. 


IV. CORRELATIONS IN FR-RMF MODELS 

In the context of QHD, protons and neutrons are 
the fundamental particles interacting each other through 
scalar and vector mesons exchange. In this framework, 
the fields a and to represent, respectively, these mesons 
and mimic the attractive and repulsive parts of the nu¬ 
clear interaction. The main representative of QHD mod¬ 
els is the Walecka one [5] , in which the only two free pa¬ 
rameters are fitted in order to reproduce the values of p 0 
and B a . However, it does not give reasonable values for 
K a (~ 500 MeV), and M* (~ 0.54M). This problem was 
circumvented by Boguta and Bodmer [13], who added to 
the Walecka model cubic and quartic self-interactions in 
the scalar field a, introducing, consequently, two more 
free parameters, which are fitted so as to fix these quan¬ 
tities. All thermodynamic quantities of this model are 


found from its Lagrangian density, given by 
£ = - M)ip + g a (ripip + ^{d^ad^a - m 2 a a 2 ) 

- j* 3 - |a 4 - gjpj^ip - ^F^ v F pv + 

- y VD'Vm'tV’ - i m 2 p ppp», (54) 

with F pv = d v tjj p - d p u v , and B pu = d v p p - d,,p v - 
g p (p p x p v ). The coupling constants are g a: g' u , g p , A 
and B. For a complete description of the model, and 
also other kind of FR-RMF ones, such as density depen¬ 
dent, crossed terms and nonlinear point-couplings, we 
address the reader to the recent study of Ref. [6] in¬ 
volving an analysis of 263 relativistic parametrizations 
under constraints related to symmetric nuclear matter, 
pure neutron matter, symmetry energy, and its deriva¬ 
tives. Here, we mainly focus in searching for correlations 
between bulk parameters for the FR-RMF parametriza¬ 
tions described by Eq. (54). 


A. Isovector sector 

Besides its complete analytical structure, other advan¬ 
tage of the NR limit of NLPC models is its usefulness 
in predictions of correlations in Boguta-Bodmer models, 
as pointed out in Ref. [12]. For instance, the L 0 x J 
correlation in Eq. (26) is showed to be linear also for 
these models under the restriction of fixed values of ef¬ 
fective mass. Like in the NR limit, different values for 
K a do not destruct the linear dependence, see figure 2b 
of Ref. [12]. Based on this correspondence, we use the 
framework of the NR limit to confirm other correlations 
in Boguta-Bodmer models. Still at the isovector sec¬ 
tor, we showed in Ref. [12] that the linear correlation 
indicated in Eq. (28) holds for the models described by 
Eq. (54), if we also fix the values of K 0 and J. Now, 
we further investigate such a correlation. In Fig. 17, 
we show K° as a function of L 0 for parametrizations 
with effective masses submitted to the FRS constraint, 
Eq. (1). According to Ref. [10], this is the range of 
to* in which Boguta-Bodmer models have to be con¬ 
strained in order to produce spin-orbit splittings in agree¬ 
ment with well established experimental values for the 
16 0, 40 Ca, and 208 Pb nuclei. In this figure, we present 
curves corresponding to the limiting values of the ranges 
25 ^ J < 35 MeV [ 6 ], and 250 < K 0 < 315 MeV [29], 

From these results, we can conclude that the relation 
K° m = PreiL 0 + Qrei also works well for the Boguta- 
Bodmer models submitted to the FRS constraint. How¬ 
ever, by comparing the angular coefficients p(p 0 ) from 
Eq. (28) and p re i, we notice that p re i slightly depends 
on K a , unlike the nonrelativistic case in which the an¬ 
gular coefficient depends only on p a . We perform the 
same analysis for the L 0 dependence on Q° yra , motivated 
by the linear correlation presented in Eq. (31). The re¬ 
sult is found in Fig. 18. As we see, there is a corre- 
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FIG. 17. (Color online) K° ytn versus L 0 for Boguta-Bodmer FIG. 19. (Color online) Symmetry energy slope as a function 
models with p 0 = 0.15 fm~ 3 and B 0 = 16 MeV. The linear of p/po for some Boguta-Bodmer parametrizations. 
fittings are indicated in both panels. 



FIG. 18. (Color online) Q° yrn versus L 0 for Boguta-Bodmer 
models with p 0 = 0.15 fm~ 3 , B 0 = 16 MeV, and 0.58 ^ m* ^ 
0.64. The linear fittings are indicated in both panels. In the 
insets, we show L 0 in the ranges of (a) 79 ^ L 0 ^ 83 MeV, 
and (b) 109 ^ L 0 ^ 113 MeV. 


lation between Q° ym and L a . However, the linear form 
Q°ym = u reiL 0 + v re i is strongly dependent on K a , unlike 
the previous case. For K a = 315 MeV, for example, we 
notice the range of effective mass that still ensure a line 
in the curve fitting is reduced from the range of Eq. (1) 
to 0.60 ^ to* ^ 0.64. The break of this linearity is better 
depicted in the insets of Fig. 18. 

One can also use the linear dependences showed in 
Figs. 17 and 18, with the latter guaranteed only for some 
values of K 0 and to*, in order to justify a possible crossing 
point in the density dependence of the symmetry energy 
slope, as we did in the case of nonrelativistic models in 
Sec. III. In fact, there is such a crossing as we can see 
in Fig. 19. From the angular coefficients of the respec¬ 
tive lines of Figs. 17a and 18a, we can solve the equation 
1 + 8.87x% + = 0 to find a crossing density at 

Pc/ Po = 0.61. This value can not be refined since a linear 
correlation is not found in the next order bulk parame¬ 


ter, namely, I° ym . as we see in Fig. 20. We see that a 



FIG. 20. (Color online) 7° m versus L 0 for Boguta-Bodmer 
models with p 0 = 0.15 fm -3 , B 0 = 16 MeV, and 0.58 ^ m* ^ 
0.64. 


correlation between L a and I° ym holds for the relativistic 
Boguta-Bodmer models, but it is not a linear one, as we 
verified in the NR limit, Eq. (37). 

Still concerning isovector bulk parameters, we point 
out to the reader a specif class of relativistic models of 
Ref. [35] with mesonic crossed interactions. Following 
notation of Ref. [ 6 ] , they are classified as type 4 models 
(<r 3 + a 4 + Uq+ cross terms models) and have the terms, 

£<rcjp = ^04 4" 2 011 

+ + gag^P^ ^2 + ^a 2 'g<j<^j 

(55) 

added to the Lagrangian density of Eq. (54). The 
parametrizations of this model presented in Ref. [35] are 
constructed in order to fix the symmetry energy not at 
the saturation density, but in a smaller value. They 
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present S(pl) = 26 MeV, with p s c /p 0 = 0.68. Therefore, 
the S(p) function presents a crossing point, as pointed 
out in Ref. [16], and as one can see in Fig. 21a. Therefore, 



FIG. 21. (Color online) (a) Symmetry energy as a function of 
p/po, and (b) L 0 as a function of J for the parametrizations 
of Ref. [35] . 


such a crossing indicates a linear behavior between L a 
and J for these specific parametrizations. This correla¬ 
tion is clearly observed in Fig. 21b. Furthermore, the 
crossing density is obtained from the angular coefficient, 
by solving the equation 1 + 9.27a;® = 0. The solution of 
this linear equation leads to pi/po = 0 . 68 , exactly the 
value verified in Fig. 21a. 

As a last remark of this subsection, we point out to 
the reader that the angular and linear coefficients of the 
L a x J correlation are not universal, as we can see by 
comparing the linear equation of Fig. 21b of the rela¬ 
tivistic FSU family, with that of Fig. 9a of the Skyrme 
parametrizations. Even among relativistic models, one 
can not reach such universality, as we can see by the com¬ 
parison of the correlation in Fig. 21b with that found in 
Ref. [36] for the relativistic NL3* and IU-FSU families. 


B. Isoscalar sector 

Motivated by the analytical structure relating K a , Q 0l 
and I a , we investigate in this section whether the lin¬ 
ear dependences presented in the NR limit, see Eqs. (43) 
and (47), also applies for Boguta-Bodmer models. Ac¬ 
cording to the NR limit case, if we keep fixed the ef¬ 
fective mass, K a linearly correlates with Q a as we show 
in Eq. (43) and Fig. 14a. For the relativistic case of 
FR-RMF models described by Eq. (54), we see that this 
condition remains, as one can see in Fig. 22a for the 
MS2 [37], NLSH [38], NL4 [39], NLRA1 [40], Q1 [41], Hy¬ 
brid [42], NL3 [38], FAMAl [43], NL-VT1 [44], NL06 [ 6 ], 
and NLS [45] parametrizations presenting m* ~ 0.6. 
Such a correlation can be used in order to justify the 
crossing in the K{p) function depicted in Fig. 22b. Pro¬ 
ceeding in that direction, we use the expansion of the 


i 1 r 

■ (a) 

/■ 



1 ' 1 T \ 

— MS2 

^ 1 1 1 
(b) / 



300 


-- NLSH 






■■■■ NL4 

/ / 





NLRA1 


# / 




— Ql 

/ *// 
• :U 

„r ; / • 


200 


-- Hybrid 

/ . */*_ 

MS2 



NL3 



NLSH 



FAMAl 

1 ■ 

" JF ♦ 

NL4 



—* NL-VT1 , 


cv° 7 A 

NLRA1 

100 


NL06 U 

/ _ 

/ < 

Q1 


NLS 


/ v 

/ ► 

Hybrid 

NL3 



II 

^ O 

s 

0.77p o ' 

■ / + 

FAMAl ' 

0 



- 

. / x 

NL-VT1 . 



bjf’ / 



NL06 



f / m* - 

0.60 

-/ 

NLS - 





L _1_1_L 




_1_i_1_i_1_ 

_1_I_ 


_i_,_i_,_ _i_i_,_i_,_i_,_i_,_ 

200 300 400 0.5 0.6 0.7 0.8 0.9 1 

K(MeV) p/p Q 


FIG. 22. (Color online) (a) Q 0 as a function of K a for some 
Boguta-Bodmer parametrizations with m* ~ 0.6. (b) Incom¬ 
pressibility as a function of p/po for the same parametriza¬ 
tions. 

energy per particle in Eq. (50) until order a; 3 to calcu¬ 
late the density dependence of the incompressibility. The 
result of this calculation is given by 

K{p) ~ (3s + 1) [K 0 + (9 K a + Q 0 )x + 6 Q a x 2 } . (56) 

Therefore, the linear dependence Q 0 = i re iK 0 + j re i 
showed in Fig. 22a can be used in Eq. (56) to furnish 

K(p) ~ (3s + 1) {Ko [l + (9 + i re i)x + 6is 2 ] 

+ (1 + 6s) jreix} , (57) 

with i re i = 6.30 and j re i = 1552 MeV. Thus, one has 
a crossing point when the quadratic equation 1 + (9 + 
VeO^c + 6* re ;s]f = 0 present solution. This is the case 
for Pc/po = 0.78. 

We remark that a crossing point in the K(p) func¬ 
tion was firstly explained from the linear correlation be¬ 
tween K a and Q a in Ref. [14]. However, the authors 
found such a crossing for some nonrelativistic Skyrme 
and Gogny parametrizations. For the relativistic mod¬ 
els analysed, they did not found linear correlations or 
crossing points. Indeed, for the nonrelativistic case, they 
found a crossing at p*/p 0 ~ 0.7, a value quite close to 
ours findings, namely, p^/po = 0.77 and p^/po = 0.79, 
for Boguta-Bodmer models and the NR limit, respec¬ 
tively, see Figs. 22b and 15. 

Unlike the linear correlation presented in the NR limit, 
the angular coefficient i re j is slightly dependent on the ef¬ 
fective mass. For the former case, we have i = i(p a ), 
see Eq. (44). In Fig. 23, we show this variation ob¬ 
serving the FRS constraint and the range 250 ^ K a ^ 
315 MeV. In particular, notice that for m* = 0.64, a 
value that ensure good values for finite nuclei spin-orbit 
splittings [10], the range of the skewness coefficient is 
given by —183 < Q 0 < 130 MeV. Such a specific con¬ 
straint for Q 0 present an overlap of ~ 36% with a re¬ 
cent range proposed for this bulk parameter in Ref. [46], 
namely, —494 ^ Q 0 ^ —10 MeV. In this study, the au¬ 
thors analysed models with crossed interactions among 






14 



FIG. 23. (Color online) Q 0 as a function of K 0 for Boguta- 
Bodmer parametrizations in which p 0 = 0.15 fm” 3 and B 0 = 
16 MeV. The linear fittings are indicated in the figure. 


the fields, i. e., models described by Eq. (54) added to 
the terms in Eq. (55). They verified that such mod¬ 
els, presenting the skewness coefficient in the range of 
—494 ^ Qo ^ —10 MeV, satisfy the suprasaturation 
constraint for the density dependence of the pressure in 
the symmetric nuclear [47], and also the neutron star 
mass constraint, given by 2.01 ± O.O4M 0 . This lat¬ 
ter is due to the recently discovered neutron star PSR 
J0348+0432 [48]. 

Finally, we verify whether the relationship between K a 
and I 0 presented in Eq. (47) is preserved in the relativistic 
case. According to the results of Fig. 24a, we see that 
the linear behavior I 0 = w re iK 0 + z re i still remains in 
the range of 250 ^ K a ^ 315 MeV, with w re i slightly 
depending on m*, like in the case of i re i. 



FIG. 24. (Color online) I 0 as a function of K 0 for Boguta- 
Bodmer parametrizations in which p 0 = 0.15 fm” 3 and B 0 = 
16 MeV. The curves are constructed in the range of (a) 250 ^ 
K 0 ^ 315 MeV, and (b) 180 ^ K 0 ^ 360 MeV. The linear 
fittings are indicated in panel (a). 

It is worth to notice that, as showed in Fig. 24b, for 
a broader range of K 0 the linear dependence is blurred, 


although a correlation between the bulk parameters K a 
and I 0 still remains. 


V. MODELS WITH MORE THAN ONE 
ISOVECTOR COUPLING CONSTANT 

In previous sections, we have analyzed under what con¬ 
ditions the linear correlations presented in the NR limit 
of point-coupling models are reproduced in the context 
of relativistic Boguta-Bodmer parametrizations. How¬ 
ever, our comparisons were restricted to relativistic and 
nonrelativistic models presenting only one isovector cou¬ 
pling constant, namely, g p and G^ v , respectively. For 
the Boguta-Bodmer model, g p is related to the interac¬ 
tion strength between the nucleon and the p meson. For 
the NR limit model, G^ v regulates the strength of the 
term that mimics the same kind of interaction (we re¬ 
mind the reader that our NR limit model is derived from 
a relativistic point-coupling model, therefore, a model 
in which there are no meson exchanges). Regarding the 
specific relationship between L a and J, we concluded in 
Ref. [12] that for the NR limit model, such correlation 
is linear whenever the isoscalar bulk parameters, namely, 
to*, p 0 , B 0 and K a , remain unchanged. We also showed 
that this same condition also ensures a linear correlation 
between L a and J for Boguta-Bodmer parametrizations. 
Furthermore, the angular coefficients of these correlations 
are the same, and the absolute values of L a are very close 
each other, as we can see in Fig. 25 for the NL3* Boguta- 
Bodmer parametrizations and their respective NR limit 
versions, namely, the NR-NL3* ones. 



FIG. 25. L 0 as a function of J for the NL3* Boguta-Bodmer 
parametrizations and their NR limit versions NR-NL3*. 
These parametrizations were constructed by fixing p 0 = 
0.15 fm -3 , B 0 = 16.31 MeV, m* = 0.594, K a = 258.25 MeV, 
and by running J. 

In order to construct this figure, we fixed the isoscalar 
parameters values of the NL3* model, and varied the J 
values. We taken such procedure for the exact relativistic 
NL3* parametrization, and for its NR limit version. For 
the latter, we have used our Eq. (26). 
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By proceeding one step further in our analysis of the 
L 0 x J correlation, we now study relativistic and non- 
relativistic models with more than one isovector param¬ 
eter in order to verify whether the dependence observed 
in Fig. 25 still applies. For the relativistic model, we 
use that described by the Lagrangian density of Eq. (54) 
added to those of Eq. (55) with a.\ = a .2 = ot\ = a' 2 = 
C = 0, i. e., we choose a model with interaction between 
the mesons w and p. Thus, we are dealing with a model 
with two isovector parameters, namely, g p and a 3 . 

To take the NR limit of this specific model, we con¬ 
struct the following point-coupling Lagrangian density, 


£ nlpc = V’Ob''^ - M)i/j - ^GU^i/j) 2 + ^G 2 s {4>ip) 2 
+ ^(#) 3 + jGWO 4 - ^Glvitp^Ti/j ) 2 
- iGV TV (^7^V0 2 (^7^) 2 , (58) 


in which the last term mimics the interaction between 
the mesons w and p. The isovector coupling constants 
of this model are G 2 V and G vTV . The symmetry energy 
and its slope for the NR limit of this model are given by, 


iS (NR) (p) — G^vP + GvtvP 3 + 


\pl 

6M*(p, 1/2) ’ 


(59) 


and 


^ (m + \gI„ + SAP +1 b/) 

+ 3G 2 v p + 9G VT vP 3 , (60) 


respectively. The new isovector coupling constant, G v tv, 
is found by imposing upon the model that the symme¬ 
try energy at pi/p 0 = r is fixed at a particular value 
Si = 5(|0i). Here, r is a value smaller than 1. Further¬ 
more, we still found G 2 V by requiring that the model 
present a particular value J for the symmetry energy at 
the saturation density. The analytical form of these con¬ 
stants as a function of the bulk parameters can be found 
in the Appendix. 


with 


b' (to*, p 0 , B 0 , K 0 , r, Si) = b(m*,p 0 , B a , K a ) 
10 E°r 2 / 3 

+ 3(r 2 - l)(3Af 2 - 19 E°M + 18 Ef) * 
x{^ (13M - 66E P °) + ^ (M - 6E°) 

- B a (9 M - 48 E°) - E° (7 M - 36 E°) 

(M - 6 E°) + Ik o (M - 3 El) 


- 6 B 0 (3 M - 13 E°) - 2 E° (5 M - 27 E°) 


|f (M - 10E f °) - ^ (M - 2E°) 


+ 3 B 0 (3 M - 10 E°) + 3 E° (M - 6 E°) 


6 


0 r 2 - 1 ) 


Si + 5£g _ r 2 /3 . + 5E^ , 2/3 _ r -l/3\ 

r 9 to* 9 


(62) 


Notice that now, a linear correlation between L 0 and J 
is established if the function b'(m*,p 0 ,B 0 ,K 0 ,r,Si) is a 
constant, i. e., if the quantities to*, p a , B 0l K 0 , r, and Si 
are kept fixed. Moreover, if we now look at the S(p) func¬ 
tion for a particular parametrization family, namely, that 
in which the set to*, p 0 , B a , K 0l r, and <Si is fixed and J 
runs a certain range, we see a crossing point, differently 
from the NR limit case presenting only one isovector cou¬ 
pling constant. We show this finding in Fig. 26 for the 
NR-NL3* family. 



FIG. 26. Density dependence of the symmetry energy for the 
NR-NL3* parametrization family with two isovector parame¬ 
ters. 


Such an analytical structure enables us to find the fol¬ 
lowing correlation between L 0 and J , 


L n = 3 


3 — r 


J + b'(m*,p 0 ,B 0 ,K 0 ,r,Si), (61) 


By looking at the finite range relativistic model with 
the ui and p mesons interaction, we verified that a linear 
correlation between L a and J also holds if we apply the 
same conditions observed in the NR limit case, i. e., fixed 
values of to*, p a , B a , K a , r, and <Si. A direct compar¬ 
ison between these two correlations, analogous to that 


1 — r 2 
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presented in Fig. 25, is displayed in Fig. 27. Here, we 



FIG. 27. L 0 as a function of J for the relativistic NL3* 
parametrization family and their NR limit versions NR-NL3*. 
These parametrizations were constructed by fixing p 0 = 
0.15 fm“ 3 , B 0 = 16.31 MeV, m* = 0.594, K a = 258.25 MeV, 
r = 0.67, <Si = 25.8 MeV, and by running J. Both models 
present two isovector coupling constants. 

restricted our analysis for J in a range of values greater 
than Si. 

From Fig. 27, we can notice that the NR limit ver¬ 
sion of NL3* parametrizations with two isovector cou¬ 
pling constants, presents a different slope for the L a x J 
linear correlation, differently from the case showed in 
Fig. 25, where we tested models with only one isovec¬ 
tor parameter. This result is in qualitative agreement 
with the findings obtained in Ref. [36], where the au¬ 
thors compared the same NL3* parametrization family 
(two isovector parameters) with a nonrelativistic Skyrme 
parametrization family named as SkNL3*. For this fam¬ 
ily, the isoscalar bulk parameters present the same values 
as in the relativistic NL3* model. The authors also im¬ 
posed that the energy per neutron predictions, at subsat¬ 
uration densities, of the SkNL3* and NL3* models were 
compatible with the band constraint depicted in Fig. 2 
of Ref. [36]. As a consequence, they found correlation 
bands (ellipses) for the L 0 and J bulk parameters, as one 
can see in Fig. 28. 

These ellipses were constructed by the authors of 
Ref. [36] from the covariance analysis method. In Fig. 28, 
we extracted such bands and constructed the linear fits. 
Notice the nonrelativistic line presenting a greater slope 
in comparison with the relativistic one, exactly the same 
qualitative behavior observed in Fig. 27, where we have 
constructed the correlations only observing the condi¬ 
tions under which they are linear ones. Furthermore, 
our ratio for the NL3* slope to the NR-NL3* one is not 
much different to the same ratio of Fig. 28, namely, 2.22 
for ours (Fig. 27), and 2.96 for Ref. [36], or Fig. 28, ob¬ 
tained from the covariance analysis method. 



FIG. 28. (Color online) Correlation bands for the L 0 and 
J bulk parameters of the NL3* parametrization family and 
SkNL3* one, extracted from Fig. 3a of Ref. [36]. The full 
lines are linear fits. 


VI. SUMMARY AND CONCLUSIONS 

In this work, we analysed the arising of correlations be¬ 
tween isovector and isoscalar bulk parameters of hadronic 
nonrelativistic, and relativistic mean-field models. In 
particular, we discussed the connection of the crossing 
point in the density dependence of a particular bulk 
quantity, with the specific linear correlation between this 
quantity with its immediately next order bulk parame¬ 
ter. In the isovector sector, for instance, if there is a 
crossing point in the density dependence of the symme¬ 
try energy, then, it can be explained by the linear cor¬ 
relation between the symmetry energy, S, and its slope, 
L = 3p(dS/dp ), both evaluated at the saturation density, 
i. e., there will be a linear correlation between J = S(p 0 ) 
and L a = L(p 0 ). In summary, the crossing points can 
be seen as a signature, or a route, in the searching of 
linear correlations among bulk parameters, as discussed 
in Sec. II. 

In the nonrelativistic framework, we presented corre¬ 
lations in some Skyrme [3] and Gogny parametrizations, 
see Figs. 2, 7, 9, 10 and 14b, as well as in parametrizations 
generated from the NR limit of NLPC models. By us¬ 
ing the analytical structure of the latter model, we could 
write its five coupling constants in terms of the bulk pa¬ 
rameters p 0 , B a , K a , m*, and J in order to investigate 
the conditions in which the linear correlations of L a with 
K° , Q° and 1° , in the isovector sector, and of K 0 
with Q 0 and I 0 , in the isoscalar one holds. For these 
parametrizations, we showed that L a linearly correlates 
with K° „,Q" and 1 ° m if we keep fixed the values of J 
and K a , see Eqs. (28), (31) and (37). Following analogous 
procedure, we found that parametrizations with fixed ef¬ 
fective mass lead to linear correlations of K a with Q a and 
I 0 , according to Eqs. (43), (47), and the respective subse¬ 
quent discussions. For some of these linear correlations, 
we discussed how they could have been found from the 
searching of crossing points in the bulk parameter as a 
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function of the density. We pointed out that the cross¬ 
ing at Pc/po = 0.47 ( Pc/Po = 0.79) exhibited in Fig. 3 
(Fig. 15) for the L (NR> (ii" (NR) ) function, for instance, is a 
signature of the linear correlations between L 0 and K° ym 
( K a and Q 0 ), at least. 

Regarding the relativistic mean-field models [6], we 
mainly studied that presenting cubic and quartic self¬ 
interaction in the scalar field a, namely, the Boguta- 
Bodmer model [13]. The reason for this choice was based 
on our previous work of Ref. [12]. In that work, we 
showed that some correlations among bulk parameters 
presented in the NR limit of NLPC models, are also 
valid for this particular relativistic model. We further 
studied the correlations for bulk parameters of isovec¬ 
tor and isoscalar sectors, mainly in the ranges of ef¬ 
fective mass, symmetry energy, and incompressibility 
given by 0.58 ^ to* ^ 0.64, 25 ^ J ^ 35 MeV, and 
250 K a ^ 315 MeV, respectively. The first range was 
proved to be experimentally consistent with finite nuclei 
spin-orbit splittings, according to Ref. [10]. The second 
is compatible with experimental values from analyses of 
different terrestrial nuclear experiments and astrophysi- 
cal observations [6, 28], and the latter was based on the 
recent reanalysis of data on isoscalar giant monopole res¬ 
onance energies [29]. 

In the isovector sector, we showed that L a also cor¬ 
relates with K° , Q° and 1° like in the NR limit. 
However, the linear behavior of L a with I° yrtl is broken 
in the range of m* and K 0 analysed, as pointed out in 
Fig. 20. For the correlation between L a and Q° ym , we 
verified that the linear behavior is blurred only at higher 
values of K 0l see Fig. 18. We still concluded that the 
linear dependence of L 0 on K° is preserved for fixed 
values of J, as in the case of the NR limit, according to 
the results presented in Fig. 17. This specific linear cor¬ 
relation was used to justify the crossing point exhibited 
in Fig. 19 for the L(p) function. 

By comparing the behavior of K 0 and Q 0 in the 
isoscalar sector, we verified that these quantities are lin¬ 
early correlated if the effective mass is kept fixed, exactly 
as deduced in the NR limit case. This correlation was 
displayed in Figs. 22a and 23. We also used the angular 
coefficient presented in the former figure in order to jus¬ 
tify the crossing point in the incompressibility function of 
the parametrizations showed in Fig. 22b. Furthermore, 
we notice that for to* = 0.64, and 250 ^ K a ^ 315 MeV, 
Qo varies in the range of —183 < Q 0 < 130 MeV, 
and present an overlap of about 36% with the range of 
—494 ^ Q 0 < —10 MeV recently proposed in Ref. [46]. 

Lastly, we verified that the linear behavior between K a 
and I 0 is still valid for the Boguta-Bodmer models with 
fixed effective mass and for the range of 250 ^ K a ^ 
315 MeV, see Fig. 24a. Nevertheless, the linearity is bro¬ 
ken for a broader range of K a but with a correlation 
I 0 = I 0 {K 0 ) still applying, see Fig. 24b. 
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Appendix A: Coupling constants of the NR limit 

In the NR limit of the NLPC models, the coupling 
constants can be written in terms of the bulk parameter, 
namely, to*, p 0 , B 0l K a and J, as 


G 


2 

TV 


-i 

po 6 M m* 


(Al) 


q M/3 Pi 

(3M 2 - 19E°M + 18E° 2 ) | m* K 

+ K 0 (. M - 3 E°) - 9 B 0 (3 M - 13 E°) 

- 3E° (5M - 27F;°) l, (A2) 


B = 


M/3 pi El 

(3 M 2 - 19 E°M + 18 E° 2 ) | to * 

- -E (M — 2 E°) + 3 B 0 (3 M - 10 E°) 

+ 3E°(M-6E°) l, 


(M - 10 E°) 
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- 1 - 2 A Po - 3 Bpl, 
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IL-^-Ap 0 -2Bpl. (A5) 
m* p 0 p 0 


Po \m 

In the case of the NR limit obtained from the La- 
grangian density of Eq. (58), the two isovector coupling 
constants are written as 

Si 5 E\ 


G v tv — \ J - 


(1 _ #/ 3 ) _ ^ 
9to* v ; 9 M 


M 1 - 


+ 2Ap 2 (l-r) + 3^(i_ r 2 ) 


#/3 


P 3 o(l-r 2 y 

(A6) 


and 


G'L = — - —Po i - G VTV Po, (A7) 


Po 6 M m* 

with r = pi/po and <Si = S(p\). 
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